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ABSTRACT 

A need for more knowledge of meteorological parameters in 
the region 30,000 to 45,000 feet has come about in recent years 
as a result of rapidly growing jet aircraft operations. This 
paper discusses one of these parameters, the height of the 
level of maximum winds, which affords considerable aid to the 
flight planner in his selection of preferred flight altitudes 
os The deviations of the level of maximum winds from 
the tropopause are discussed and graphically shown, indicating 
that the latter is not a satisfactory first approximation to 
height of the former. This fact gives rise to the development 
of a series of simple and multiple regression equations for 
determining the height of ae lever of naximum winds ag a 
function of both the tropopause and Eve 300-millibar heights. 

Professor Robert J. Renard of the U. S. Naval Post- 
graduate School has been the guiding force behind this investi- 
gation. For his time and effort, the writer is deeply 


indebted. 
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L. inm@roduction. 

With the advent of "yet%aircraf® in” recent years, the Pore- 
caster‘s spatial envelope of responsibility has almost doubled. 
It has become necessary to analyze the meteorological parameters 
up to 100 mb on a daily basis. Knowledge of wind, cloud, tem- 
perature, and pressure fields at these levels are necessitated 
by the demands of aircraft flight planning. The United States 
Weather Bureau [97\ has suggested that the forecast accuracy of 
the temperature and wind at jet flight levels be + 3°C and 
+ 20 knots respectively. 

Not the least important of the parameters needed for econom- 
ical operation of aircraft over long distances are the wind 
speed and direction, including their variation both in time and 
space. Over the continental United States and parts of Canada, 
a United States Weather Bureau analysis of the layer of maximum 
wind (hereafter abbreviated LMW and used synonymously with level 
of maximum wind) is available on a 1/40,000,000 polar stereo- 
graphic facsimile chart every 12 hours. The height of the LMW 
typically varies from about 24,000 feet in the northern sec- 
tions “to about 44,000 feet in the south. This analysis, along 
with the contour analyses of the standard isobaric surfaces, 
reasonably describes the horizontal wind field and its deriva- 
tives throughout the high troposphere and lower synoptic strato- 


sphere.* However, the extension of the LMW analysis, particularly 


* a term defined by H. A. Panofsky fol » Which includes the 
region from the tropopause to about 10 millibars. 
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to ocean areas, can be accomplished only with great difficulty 
owing to the paucity of wind data in the stratosphere. Thus, 
some means of determining with reasonable accuracy, the LMW data 
from the more easily obtained variables at lower levels would 
allow extension of the high-level wind analysis to ocean areas. 
One phase of the problem involves determination of the height 

of the LMW. It is this parameter which is under statistical 


investigation here. 





2. Background, 

Much has been written concerning the environment of the jet 
aircraft since the inauguration of high altitude, lony range 
flights within the past ten years. Relative to the specific 
problem at hand, Reiter (5 \ discusses many aspects of the iayer 
of maximum winds, including reliability of upper-wind data, 
analysis of the horizontal and vertical structure of the LMW, 
and forecasting the LMW, with special application to jet flight 
operations. His general conclusions testify to the need for not 
only more accuracy in upper -wind measuring devices, but to theix 
proper representation and interpretation as well. 

Since there are rather large inaccuracies in the upper- 
wind data and because of the scarcity of data in certain areas, 
it becomes desirable to determine the LMW from more commonly avail- 
able data, The United States Weather Bureau assumes that, for 
purposes of analysis, the height of the LMW is at, or very nearly 
at, the tropopause in middle and northern latitudes, The tropo- 
pause, if considered a first-order discontinuity in the tempera- 
ture field, is logically where one would expect to find the 
maximum wind. Unfortunately, as will be shown later, there are 
many occasions in which the LMW deviates significantiy from the 
tropopause height. Johannessen (2a) discusses a series of ob- 
servations taken at Larkhill, England, which indicates that 
the LMW occurs from 2000 to 4000 feet below the tropopause. 
However, extension to other latitudes is not plausible. 


The investigative results whose discussion follows is a 
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consequence of the actual deviations of the height of the LMW 
from that of the tropopause over the United States during the 
months of November 1959 and January 1960, which include a 


great variety of weather types. 





3. Deviations of the Height of the LMW from the Height of the 
Tropopause 


The tropopause definition as used by the United States Weather 
Bureau {11} is, using only the lapse-rate criteria, that point in 
the sounding at which the lapse rate decreases to 2°C/km and then 
averages less than this in the 2-km layer immediately above. 

Whether nationally or internationally defined, the tropopause 
level is a controversial subject. Stinson [3] in recently proposing 
a new tropopause identification and classification system deals 
with a aad of soundings which bring out the difficulties in. 
selecting a certain point fulfilling the above lapse-rate condition. 
ogee less of inaccuracies or uncertainties of ae 
location, the heights as given on the United States Weather Bureau 
00Z and 122 facsimile charts of tropopause heights (27) were used, 
These charts were analyzed by the author for use in this investiga- 
tion. A typical example of the analysis is shown as part of 
Figure l. A striking feature is the band of strong height eradient 
oriented generally east-west from southern California through 
Texas to Florida. This represents, for the purpose of this study, 
a continuous zone of transition between the tropical and middle- 
latitude tropopauses. Undulations of this band are always present 
in varying degrees of amplitude and are associated with the major 
weather systems. For example, the passage of a low or trough at 
300 mb is synonymous with a trough in the tropopause heights. 
Another zone of large height gradient appears in northwest United 


States and is probably associated with the transition from the 
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middle to the polar tropopause, A hemispheric analysis o1 
this same phenomenon [3] shown in Figure 2, clearly exhibits 
the nature of the two transition zones justi mentioned. 

As an approach to defining statistically some relations 
between the height of the LMW (Zz)c)) and the tropopause, the 
United States Weather Bureau analysis of the height of the LMW 
(see Figure 1) is compared to the locally produced trepopause 
Be iysis by a graphical subtraction. of the contour heights of 


the two surfaces.’ This also is shown in Figure 1, The generai 


4 


ed 


area of approximate cotncidence of the two surfaces extends 
from. extreme northwest United States toward the southeast 
chromen Missouri and from extreme northeast United States 
toward the Pitiwes: cee Tennessee. In general, this area 
of coincidence, defined as that aa ten in which 2 Zig EO, 
follows closely the 36, 000 to 40,000-f£t channel in the tropo- 


pause field, In areas of tropopause height less than 36,000 


faer , Zanwe? Er while tropopause heivhts ‘reater. than 


= 
7 


40,000 feet are associated with By <2e '.. Further, since 
Be tropepaumt is closely associated vith-major weather systems, 
it 1s reasonable to.suspect that the’ ‘troughs , ridges, and ae 
Senbero an Figure 1 will have continuity from day to day. That 
this is true is shown in Figure 3, which is a graphical sub- 
traction for the situation as it existed 36 hours previous 


to that in Figure 1, The movement of the minus eight isoline: 


over Montana corresponds closely to the movement of a 


27,600-ft low at the 300-mb surface situated about 300 miles 





Fig. 2 


Ni 
Cie 
INS 


yi 


‘ i 

ERR O 
ms “Ue, 
o; 


7 Ft it ters 





Tropopause Height Contours Jan. 1, 1956. 
(After Defant and Taba.) 





to the north. The plus 16 area moved east and increased to a 
plus 20 over Texas with another plus 16 coming in over central 
California on Nov 17. These large plus areas appear to follow 
the movement of a 300-mb ridge located 400 miles to the west. 
Thus, figures 1 and 3 show that an assumption of ay equalling 
Zrmw can be made only over a ey ae of the continental 
United States. For the two months rae the tropopause 
was found to seit from 12,000 feet below the LMW in northern 
United States to as high as 24,000 feet above the LMM in the 


% 


south. 
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4, Statistical Relationship Between the Height of the LMY and 
other Parameters 


A study of Figure 1 suggests a definite relationship between 
(Zp-Zyyqy) and Zp. A better picture of this relationship is 
shown in Figure 4 as a scatter diagram of Zp versus (Zy7~-Zy)q)), 
the points being randomly selected from November 1959 and January 
1960. A best-fit line may be drawn through these points using 
the least squares method but because of the decrease in slope for 
the envelope of data below Zp = 36,000 feet, the standard error 
would undoubtedly be excessive. Thus, two straight lines may be 
drawn, one for Z, greater than 36,000 feet and the other for zy 
less than 36,000 feet. The discontinuity which occurs-between 
the two regression lines at Zy = 36,000 feet can easily be smoothed 
in to provide one continuous curve. Before proceeding on computa- 
tional work, approximately 1000 actual data points, distributed 
quite evenly over the two-month period, were selected, and data 
for Zp versus Zyyy recorded. About half of these randomly- 
selected points were plotted on the scatter diagram of Figure 5. 
This diagram shows that a two-dimensional linear relationship 
will give best results for values of aq ranging from 30,000 fect 
to 44,000 feet. Above and below this range the slope of the | 
data envelope changes and the scatter Mceteweee 

With a view orate simple utilization of all data points, 
it was felt that the best course of action was to work first with 
the data as originally divided, that is, for tropopausc heights 


above and below 36,000 feet. In order to use the 1000 points 
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versus yy, 
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Figure 5. 





taken over the period, it became necessary to employ the 
National Cash Register 102A electronic digital computer at the 
United States Naval Postgraduate School. The computer was 
programmed to accept a maximum sample size of 256. Outputs of 
the mean, standard deviation, variance, and linear correlations 
between the two sets of data were recorded. The correlation 


coefficient (r) used is defined in Hoel [1] as 


‘a a (X: ~ x Ome) 
é=i nN Sx os 


where sx and sy are. the standard deviations of the two varia- 
bles involved. The first correlation attempted related the 
independent variable, Zp less than 36,000 feet, to the dependent 
variable, (Zp > Zjyw), henceforth referred to as AZ. A 

sample size of 256 randomly selected points was used., This 
group of data was part of the 1000 points previously mentioned, 
A correlation coefficient of 0.66 resulted. The form of the 


regression equation, utilizing r is , 


wiIere = Syz 


od es = 2 — 
Ss 08°" eters aces 


The resulting equation is 


aZ= "I56@ Cr -IES CZ2, B/W 00's St) 


Solving for Ziywy gives 

Zen = » 494 Zr HIGH Ct) 
The standard error is 2490 feet and can be interpreted as the 
Zimw distance above and below the regression line which includes 


about 68% of the data points. The standard error defined by 
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ranotisky [a] S civen by Chevexpression 
5 ‘i 
S..E, = (/- (ame) 

During the course of the graphical subtraction mentioned 
earlier, it was noted that there exists a marked resemblance 
between the contours of the LMW and the 300-mb chart. Thus, for 
each of the 1000 points ovér the two months‘ period, a 300-mb 
height was determined and plotted against Zjmy on the scatter 
diagram shown in Figure 6, The relation here is about the 
same as that for Zy7 versus OZ. However, the 300-mb versus 
Zp plot does not show the radical change in slope associated 
with the scatter diagram of Zy plotted against Zyjy). Using 
only those points associated with Zp less than 36,000 feet, 
the correlation coefficient, ryyyw 399, is 0.67. The regression 


equation, computed as before, is 


Lamu a: Ae iil, La2c0- F644 CLrene , A 30s 7 100% ff) (2) 


with a standard error of 2230 feet. The graphs of equations 
(1) and (2) are shown in Figure /7. 

Tf the two independent variables Zp and 2399 do not account 
for the same fraction of the variance of Z;yy, there will be 
some improvement in accuracy by writing a multiple regression 
equation of the form 

a A+ be, & 6 Fag, (3) 
where a, b, and c are constants to be determined from data in 
the sample. This best-fit plane may be determined by the 


least squares method discussed in Panofsky (47 by subtracting 
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Z'tmwo the observed value of the height of the LMW, from each 
side of the above equation, and squaring as 

(2,65 - Zi.) = (a PbS, + © S3ch -— Came ¥ 
Partially differentiating the equation with respect to a, b, and 
c successively, gives three equations with a, b, and c as the 
three unknowns, All other parameters in the equation are known 
or can be calculated from the previously mentioned results 
given by the electronic computer. The simultaneous solution of 
these three equations in this particular case is a very arduous 
task requiring that an excessive number of digits be carried 
along in order to arrive at a reasonably accurate solution, 
Since the number of digits required exceeds the keyboard capacity 
of the ordinary Marchant desk calculator, the symbolic solution 


for constants b and. as given by Panofsky fa] is more practical 


and equally as valid. These equations are as follows: 


~2 
bss Simes | lentay re — [rpepestpin (55 -s00.| GB, 
— = ace 2 
S> / F400 
ay - ; 
Cz ele | Ww. jceoO — few, ST, a (5) 
Zoe aie 2 
85, 


The constant a may then be found from the following expression: 


OMG 2O 


Thus, for a solution, the simple linear correlation coefficients 
between all the variables taken two at a time, the standard 
deviation, and the means of all the variables, must be known. 


This method makes maximum use of the work previously performed 
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ae 
— - ‘jun - 
a 4 ——=_ — ‘ 
or a> i 
AL, ~ -_ =~ he 
= 2 _ . 


=) oes 





by the computer and eliminates many hours of error-prone work 
on the desk calculator. The computer results used in calculating 


this multiple regression equation are shown in Tables 1 and 2. 





Table 1. Correlation Coefficients for z= 56,000! 


a 





Standard 
Mean Deviation 
AZ -1) 205 3,313 
Za MW 32 ,338 3,009 
2309 29 ,363 758 
Zep 31,064 3,932 


Table 2. Means and Standard Deviations of 
all variables for 2, < 34,000! (IN F7) 


Thus the regression equation (3), applicable for Zp less 


than 36,000 feet, becomes 


Lun 128 -20473 pt, U7 Zy,, Ge fr, B3eq WI 100s $L) (b) 


Here the multiple correlation coefficient, K eee 


ST 


Bc. 


is 0.713 and the standard error is 2110 feet. 


Using 4Z as the dependent variable and solving the 


resulting equation for Zywyw gives 


Ferris = SLE +. 500 2- $2.57 230-55 (ALL rtAKRtAaseles /oe's FT) (7) 
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which has a Rae: Z,,2,., of 0.760 but shows no improvement 
in the standard error because of the larger variance of the 
dependent variable Zee vonbeeedl for equations (6) and (7) 
are shommath Figures 8 and 9, 


The multiple correlation coefficient via Panofsky (4 ] 


is the square root of the following equation: 


oe _@ a 
: oe ob 
en bette TN et Joo — 2 wise T ts, Avo Vz, 390 
a i learn 
2. 2 


The quantity R* gives that fraction of the total variance of 


the dependent variable that is accounted for by the two 
independent variables, 
The standard error for a three-variable linear regression 

equation was calculated by a method given by Snedecor lo] 
and involves the use of the following equation: 

ne | UaBI DA (Been = Fons) ) 2 

7 —_ 

where n is the sample size and n-3 represents the number of 


degrees of freedom. If the sample size is large, the quantity 


: = 2 = a 

EZ imu = Leite) may be accurately represented by SCZenet~ Zs wiv) 
ees n 

the variance of the dependent variable. The equation reduces 


simply to 


_— 
S£ Sis fie : 


The quantity(1-R2) gives that fraction of the variance that 
is not accounted for by the two independent variables and is 
commonly referred to as the coefficient of alienation. 


The section of the data for which 2, is greater than 


36,000 feet has been treated in exactly the same manner as 
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heretofore described for Zqy less than 36,000 feet. Tables 3 


and 4 summarize the electronic computer results for a sample 


srzevot 256. 





Table 3. Correlation Coefficients for 2.736, coo! 


Standard 

Mean Deviation 
Aa a. 2095 53 Ol 
ZaMW 38,504 3,199 
2300 3Op).1 600 
dp G2 771 >, 60D 


Table 4&4, Means and Standard Deviations of 
all Variables for 27.7 36, coo” 


The equations ca. culated from the data of Tables 3 and 4 


and their standard errors are as follows: 


Zicrw = (G4 Sr + 3028 SE, = 2e271O 2a (8) 
Liaw SoS aaron OCT oe S/F = 2FOOD Fr- (9) 
Popa = VES ~0 BB FOIA S5n AE, = SFIS (jo) 
Sunny & AGE te UO Cp F0GB 2 xgQ Fo a Fao Li 
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Equations (8) and (11) were obtained using AZ as the dependent 
variable. The graphs of equations (8) and (9) have been drawn 
in Figure 10. Nomograms have been prepared for the multiple 
regression equations and are displayed in Figures 11 and 12, 
One additional effort was made to obtain regression equa- 
tions with smaller standard errors, This was attempted, using 
the same independent variables, by choosing data over ranges 
of values of the independent variables which have the most 
apparent linear relationship with the dependent variable. The 
points randomly selected for this correlation were taken from 
the original 1000 points and were those for which 
44 c00'7 2, P28 000! 3 FA < ID Goo! . Figures 5 and 6 
show these ranges of values to be a reasonable selection. The 


computer results using one sample of size 256 are summarized 


in Tables 5 and 6 below: 





Table 5. Correlation Coefficients for 44,000> 


Zr 7 28,0007" 3 Bae, = 50, Gc’ 
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Figure ll. 
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Standard 


Mean Deviation 
dq 36,518 3,809 
Zimy 39,169 3,302 
2300 (495773 594 


Table 6. Means and Standard Deviations for 
44,000 > 2,7 28 900" 
and 4,..< 302 600% 


rhe regression equations and their standard errors are as 


follows 
Regression Equation 
4omw = . 6/4 Zr +128 S.E,=2450 Fr (/2) 
ergy = 9.93 Fay..~ 8/8 S.E&, = 2420 -7- C13) 
ae = 2ES 4. C4/ a woy Sanu SEL = 2SPO ( (4) 


The graphs of equations (12) and (13) are drawn in Figure 13. 
A nomegre-n for equation (14) is shown in Figure 14. The value 
of the iundependent variable of equation (14) is limited by the 
fact that their coefficients are small, thus giving an almost 
herizonial plane. Equations (4) and (5) show that the values 
of b and c approach zero as the three linear correlation co- 
efficients become equal. Table (5) shows the near equality of 
the values of r resulting from this particular selection of 


data, 
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>, Conclusions and Recommendations 

The assurption that ay = Zsa is considered acceptable 
only where Z. lies between 36,000 and 40,000 feet. Above and 
below this range, the error associated with this assumption 
increases rapidly. Within the continental' United States, it is 
not uncommon for the LMW to vary from 20,000 feet below the 
tropepause in the southern states to 12,000 feet above the 
tropopause in the north. 

Ie is concluded, then, that the multiple regression equa- 
tions (7) and (11), which together are applicable for all 
_— of a will give the best first approximation to ZLMWe 
There wilt undoubtedly be weather 9 y sirens which will be 
assocxated with large deviations of Zyyy from the plane of 
enguineina the conditions under which these deviations will 
occur, however, has not been investigated in this paper. 

the regression plane associated with equation (14) is 
very nearly horizontal and as a result, the two independent 
variables hawe small coefficients and — very little effect 
on the value of the dependent variable, Zw’ Further investi- 
gation with other iaeeeenaene variables in the range of 
applicability of this equation was prevented by insufficient 
time, 

it was originaily intended to evaluate the regression 
eguarions using independent data from the month of December 


1959. However, because of the availability of an electronic 


computer at the United States Naval Postgraduate School, it 


ek: 





was felt that time would be utilized more advantageously by con- 
tinuing to search for the most significant correlations among 
the variables under consideration. In addition, it is believed 
that the variety of weather systems occurring during the ones 
of November and January, make these regression equations suf- 
ficiently general so that they may be used with confidence for 
independent samples of data. 

Were this investigation to continue, consideration of 
independent variables other than 27 and Z399 would certainly be 
warranted. Since the winds are dynamically associated with the 
thermal structure of the atmosphere, the temperature at the 
300-mb surface may well show better association with Zjyy than 
the height of the 300-mb surface. Some cons has been done in 
attempting correlations with 1000-500 mb thickness values [7] : 
It would be logical to expect that thickness values for layers 
above 500 mb would give results as good or better. 

The further consideration of Zyyy, 2p and 2399 OF any 
other associated variables should include serial correlations 


and graphical solutions of their non-linear relationships. 
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